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The invention concerns a method and a system for three-dimensional reconstruction of an image 
representing the surface contours of at least an object (1), from at least a two-dimensional view of said 
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@) RECONSTITUTION DE SURFACES EN TROIS DIMENSIONS PAR UTILISATION DE MODELES STATISTIQUES. 

@) L'invention conceme un procede et un systeme de re- 
constitution d'une image en trois dimensions repr^sentant 
les contours surfaciques d'au moins un objet (1), h partir 
cfau moins une vue en deux dimensions de cet objet prise 
aux rayons X, dans lesquels on determiner la position de la 
source (7) de prise de vues dans un referentiel de referen- 
ce, on selectionne un modele pr6defini constituant une for- 
me moyenne de Pobjet et, de facon iterative jusqu'& ce que 
les contours du modele soient tels que les 6carts entre des 
rayons de retroprojection des contours de Pimage en deux 
dimensions depuis la source et la surface du modele soient 
minimaux, on selectionne une orientation et une position du 
modele dans le referentiel de reference, puis on selection- 
ner une deformation du modele pour modifier ses contours 
en trois dimensions. 
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RECONSTITUTION DE SURFACES EN TROIS DIMENSIONS PAR UTILISATION DE 

MODELES STATISTIQUES 

La prgsente invention concerne la reconstitution d 1 ima- 
ges en trois dimensions et, plus particuliSrement, la reconstitu- 
tion des contours surfaciques d'une image k partir de vues, m§me 
incorrpletes, en deux dimensions. 
5 un exenple d » application de la prisente invention est 

la reconstitution d ! images repr^sentant des elements osseux a 
partir d f images en deux dimensions prises, par exenple, aux 
rayons X. De telles images peuvent servir, par exenple, a simuler 
une intervention chirurgicale en permettant au praticien de pr€- 
10 visualiser les positions respectives des os au niveau d'une arti- 
culation, par exenple, en cas de plastie ligamentaire ou de mise 
en place de prothSse. 

Aujourd'hui, pour permettre au praticien de visualiser 
une image en trois dimensions a partir d» images en deux dimen- 
15 sions, on est contraint d'utiliser un scanner a rayons X effec- 
tuant des tomographies. A partir de ces tonographies, on sait 
reconstituer une image en trois dimensions. Soit il s'agit d'une 
technologie dite de tomographies en trois dimensions dans 
laquelle un grand nombre de vues {de l'ordre de 200) sont prises 
2 0 au moyen d'un scanner a rayons X selon diff ^rentes positions, et 
les contours de l'os sont alors reconstitues par traitement 
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d' image a partir de ces vues au scanner. Soit il s'agit d'une 
technologie dite de tomographies en deux dimensions dans laquelle 
un grand nombre de coupes sont prises au raoyen d'un scanner a 
rayons X transversalement & l'os. On peut alors reconstituer 
5 l 1 allure et la structure de l'os. 

La technique du scanner donne de bans resultats mais 
est d'une mise en oeuvre lourde et couteuse. En effet, l'utilisa- 
tion d'un scanner permet d'dbtenir un ensemble d' images en deux 
dimensions fournissant non seulement des informations sur le 

10 contour mais Sgalement sur l'interieur de l'os. Or, dans de nom- 
breuses applications/ seule la connaissance du contour surfacique 
de l'os ou de l'objet est n£cessaire. 

Un autre exetrple d f application de la pr6sente invention 
est la recanstitution d f ossements incotrplets , par exeraple, en 

15 archeologie. Ufre reconstitution d r images en trois dimensions peut 
permettre de retrouver de fa<?on quasi-parf aite la forme origi- 
nelle de l'os mime si celui-ci est decouvert de fagon incomplete. 
Dans une telle application, des problemes similaires & ceux expo- 
ses ci-dessus en relation avec la simulation d 1 interventions chi- 

2 0 rurgicales se posent. En particulier, il est souvent utile de 

connaitre la forme de l'os sans avoir a se preoccuper de sa 
structure interne. 

Un autre inconvenient des techniques connues est 
qu'elles imposent une dose de radiations innportante pour le 
25 patient ce qui n'est pas souhaitable. Si cet inconvenient est 
moins sensible dans le domaine de 1 1 archeologie ou c'est le cout 
qui predomine, il est particulierement genant dans la simulation 
d 1 interventions chi rurgicales. 

La presente invention s f applique plus particulierement 

3 0 i. la reconstitution d ! images se rapportant & des objets identi- 

fies, c'est-a-dire, dont on connait £ 1' avarice la forme generale . 
par exenple, pour un os, il faut au prealable decider de quel os 
il s'agit. 

La prisente invention vise a proposer un nouveau pro- 
3 5 c£d§ de reconstitution d' images en trois dimensions qui pallie 
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les inconvenient s des techniques connues. L» invention vise # en 
particulier, a proposer une solution qui ne necessite pas l'em- 
ploi coQteux d'un scanner a rayons X. 

La presente invention vise egalement a proposer une 
5 solution qui soit compatible avec une exposition tninimale aux 
rayons X ou equivalents. 

L' invention vise en outre at minimi ser le nombre de vues 
en deux dimensions necessaires pour reconstituer 1" image en trois 
dimensions . 

10 Pour atteindre ces bbjets, la presente invention pr£- 

voit un procede de reconstitution d'une image en trois dimensions 
repr^sentant les contours surfaciques d'au moins un objet, a par- 
tir d'au moins une vue en deux dimensions de cet dbjet prise aux 
rayons X, caracterise en ce qu'il consiste a : 

15 determiner la position de la source de prise de vues 

dans un referentiel de reference / 

seiectionner un modele pr^defini constituant une forme 
moyenne de l'objet / et 

de fagon iterative jusqu'a ce que les contours du 

2 0 moddle soient tels que les hearts, entre des rayons de retropro- 
jection des contours de 1' image en deux dimensions depuis la 
source et la surface du module, soient minimaux : 

seiectionner une orientation et une position du module 
dans le referentiel de reference, puis 

2 5 seiectionner une deformation du module pour modifier 

ses contours en trois dimensions. 

Selon un mode de realisation de la presente invention, 
le module est obtenu a partir d'une population d» bbjets pour 
laquelle on recherche la correspondance statistique commune a 

3 0 tous les objets pour determiner une forme moyenne et les deforma- 

tions principales par rapport a cette forme moyenne, de fagon a 
disposer d'au moins un modele statistique. 

Selon un mode de realisation de la presente invention, 
les etapes de selection iteratives consistent a faire subir au 
35 modele statistique, successivement, une transformation rigide 
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modifiant sa position et/ou son orientation et une deformation 
non-rigide modifiant ses contours surf aciques . 

Selon un mode de realisation de la presente invention, 
les contours de 1' image en deux dimensions sont dbtenus automat i- 
5 quement en projetant le modele dans le plan de 1' image en deux 
dimensions, et en deformant les contours projetes de fagon & les 
faire coincider avec les points de fort gradient en niveaux de 
gris de l 1 image en deux dimensions. 

Selon un mode de realisation de la presente invention, 
10 la determination automatique des contours de l 1 image en deux 
dimensions est effectu§e de fagon iterative, chaque iteration 
§tant intercalee entre deux iterations successives des etapes de 
selection. 

Selon un mode de realisation de la presente invention, 
15 on determine, dans le referent iel de reference, des coordonnees 
en trois dimensions de points de l'objet, de fa<?on & disposer de 
points de reference suppiementaires pour les etapes de selection 
iteratives de position, orientation, et deformation. 

Selon un mode de realisation de la presente invention, 
2 0 on utilise plusieurs images en deux dimensions pour lesquelles 
les positions respectives de la source de prise de vues sont tou- 
tes determinees dans le referentiel de reference, et on effectue 
les etapes de selections iteratives en tenant compte des rayons 
de retroprojection des contours de toutes les images en deux 
25 dimensions. 

Selon un mode de realisation de la presente invention, 
le nombre d 1 images utilise est fonction de la precision souiiai- 
tee, 

Selon un mode de realisation de la presente invention, 
30 la surface du modele est constituee d ? elements de triangle, les- 
dits ecarts etant mesures par rapport a des points de certaines 
aretes constituant des generateurs du contour en trois dimen- 
sions . 

Selon un mode de realisation de la presente invention, 
35 le procede est applique £ la reconstitution des contours surfaci- 
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ques de plusieurs objets lies entre eux par des relations de 
transformation rigide et/ou elastique. 

Selon un mode de realisation de la presente invention, 
le proc6d§ est applique & la reconstitution d' images d'os. 
5 L» invention concerne egalement un systeme de traitement 

d' images, ccmportant des moyens pour la mise en oeuvre du proc6de 
de reconstitution d' images en trois dimensions. 

Ces objets, caracteristiques et avantages, ainsi que 
d'autres de la presente invention seront exposes en detail dans 
10 la description suivante de modes de realisation particuliers 
faite 4 titre non-limitatif en relation avec les figures jointes 
parmi lesquelles : 

la figure 1 reprSsente un mode de realisation d f un sys- 
tdme de re const itution d' images en trois dimensions & partir d ! un 
15 dispositif de prise de vues en deux dimensions selon la presente 
invention ; 

la figure 2A illustre la prise d'une vue en deux dimen- 
sions pour la mise en oeuvre du proc§d6 de reconstitution d* ima- 
ges en trois dimensions selon la presente invention ; 
20 la figure 2B illustre la prise de trois vues en deux 

dimensions pour la mise en oeuvre du procede de reconstitution 
selon la presente invention ; 

les figures 3A, 3B et 3C illustrent trois etapes carac- 
teristiques du procede selon un mode de realisation de la pr6- 
25 sente invention ; 

la figure 4 illustre un mode de determination de g6n£- 
rateurs de contour selon la presente invention ; et 

la figure 5 illustre un mode de determination prefgrS 
selon 1 ! invention de l'gcart entre un module en trois dimensions 
3 0 et une vue en deux dimensions. 

Les memes elements ont ete d^signes par les memes refe- 
rences aux differentes figures qui ne sont pas a l'echelle. Pour 
des raisons de clarte, seuls les elements qui sont necessaires & 
la comprehension de 1 ' invention ont ete representes aux figures 
35 et seront decrits par la suite. En particulier, les moyens de 
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traitement d f images comprenant les calculateurs informatiques 
n f ont pas ete representees et ne seront pas decrits de fa<?on 
detaill£e, car leur realisation est a la portee de l'hcnme du 
m£tier & partir des indications fonctionnelles donn£es ci-apres. 
5 Une caracteristique de la presente invention est de 

rechercher la position et 1 ! orientation d'un dbjet a reconsti- 
tuer, et dont on connait au moins une vue en deux dimensions, a 
partir d f une base de donnees contenant des modeles de cet dbjet. 
Quand plusieurs vues en deux dimensions sont utilisees, ces vues 

10 sont toutes referencees dans un meme r6f£rentiel. Ainsi, 1 ' inven- 
tion concerne la reconstitution d f une image reprgsentant un objet 
dej£ identifi§ et dont on peut disposer de modules de taille 
et/ou de forme diff£rentes. 

Une caracteristique d'un mode de realisation pr£fer£ de 

15 la presente invention est d ! utiliser au moins un modele statisti- 
que deforroable, etabli a partir de la base de donnees, pour 
reconstituer la forme en trois dimensions de 1' dbjet, Ainsi, 
1 1 invention pr6voit de definir, avant toute reconstitution, une 
base de donnees contenant des modules en trois dimensions de 

2 0 1' objet a reconstituer ou, de preference, un ou plusieurs modeles 
statistiques deformables £ partir de cette base de donnees. 

La figure 1 repr^sente une vue schgmatique d'un systeme 
de reconstitution d" images en trois dimensions selon un mode de 
realisation de la presente invention. L 1 exemple de la figure 1 

2 5 concerne la reconstitution de l f image d'un os 1 £ partir de 

radiographics en deux dimensions. En figure l, on a repr£sente 
schematiquement et partiellement le corps d'un patient p dont la 
jambe contient l'os 1 que l'an souhaite visualiser. La jambe du 
patient p (done l'os 1) est, par exemple, pos£ sur une table 2. 

3 0 Une potence 3 porte un dispositif 4 de localisation en trois 

dimensions qui peut etre optique, magnet ique, mecanique ou ultra - 
sonore et qui repdre la position et 1' orientation de multiples 
capteurs-6metteurs const itu6s, par exemple, de diodes infra- 
rouges, de reflecteurs, d'emetteurs magnetiques, a ultrasons, 
35 etc. Un tel localisateur est parfaitement connu et ne sera pas 
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detaillS plus avant . On notera simplement que, pour garantir une 
localisation correcte du patient (de l f os 1), celui-ci est g6n6- 
ralement Sgalement 6quip6 d»un capteur- §met teur 18 detectable par 
le localisateur 4. En effet, selon 1' invention, toutes les vues 
5 doivent pouvoir etre exploitSes dans un meme rSferentiel (de 
reference), qui est associS & l'dbjet que l'cn reconstitue. 

Plusieurs capteurs-6metteurs peuvent §tre fix6s sur le 
syst£me de radiologic, proche de la source 7 ( capteur - emet t eur 
19) ou proche d'un d§tecteur d f image 9 ( capteur- Smetteur 5), de 
10 fagon a reperer la position du systime de radiologie par rapport 
au referentiel de r6f£rence du capteur- £mett eur 18. 

Dans certains cas, les capteurs-emetteurs 5 et 19 sant 
dif ficilement repSrables en raison de leur eloignement ou de la 
presence d'bbjets parasites dans le champ de raesure du localisa- 
15 teur 4. Dans ce cas, un capteur- 6met teur 20 peut €tre installs 
soar la table 2 dans le champ de me sure du localisateur 4. On 
positionne alors le syst£me radiologique dans de bonnes condi- 
tions de tnesure et on repSre les positions des capteur-Smetteurs 
5 et 19 par rapport au capteur-emetteur 20 une fois pour toutes 
20 (on ne renouvelle cette Stape que si on bouge le systSme de 
radiologie dans son ensemble) . 

Le systeme radiologique etant £quip6 de codeurs angu- 
laires sur ses axes, comme cela sera detains plus loin, les 
changements de position relative du systeme radiologique sont 
25 mesur£s 3. l'aide de ces codeurs et peuvent ainsi §tre reportSs 
dans le r6f€rentiel du capteur-Smetteur 20. Pour chaque prise de 
vue, on me sure la relation gecmetrique entre les capteurs- 
emetteurs 20 et 18 et, par ce biais, on reporte la position du 
systdme radiologique dans le referentiel de reference du capteur- 
3 0 emet teur 18. Ainsi, toutes les radiographics sant calibrSes dans 
le referentiel du capteur-emetteur 18 qui peut etre mobile. 

Le dispositif 6 est, par exenple, constitu§ d'une 
source 7 Sl rayons X portee par une premigre extr<§mite d'un bras 8 
en demi-cercle dont I 1 autre extr£mit6 est destinSe Si recevoir le 
35 film 9 d 1 impression de la radiographie , ou un capteur Slectroni- 
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que equivalent tel qu'un amplif icateur de brillance, ou un detec- 
teur plat au silicium amorphe. La relation entre le bras 8 et la 
table 2 est telle que cette dernidre se trouve entre la source 7 
et le capteur 9. Le bras 8 est monte & rotation autour d f un axe 
5 10, motorise ou a deplacement manuel, et porte par une potence 
11. L 1 ensemble de prise de vues en deux dimensions peut ainsi 
tourner autour de l'os 1 pour effectuer le norribre d§sir£ de 
radiographics de celui-ci. Le cas echeant, le localisateur opti- 
que 4 peut §tre associS £ un dispositif de codage angulaire de la 

10 position de I 1 axe 10. 

L f ensemble est pilote par un systeme informatique, par 
exemple, un ordinateur 12 associe, de preference, a un ecran de 
visualisation 13. A la figure l, on a symbolist par des liaisons 
unifilaires 14, 15, 16 et 17, les bus informatiques d l 6change de 

15 signaux §lectriques de commande et de donnees entre 1 1 ordinateur 
12 et, respectivement, le localisateur 4, la source 7, le moteur 
ou le codeur qptionnel de l'axe 10 et le capteur 9. 

Les figures 2A et 2B illustrent deux modes de prises de 
vues selon deux modes de realisation de la prgsente invention. 

20 Aux figures 2A et 2B, les vues en deux dimensions ont ete repre- 
sentees par les plans respectifs dans lesquelles elles sont pri- 
ses, c'est-a-dire par 1" allure de la surface du capteur 9 (figure 
1) lors des prises. Cela correspond aux images radiographiques en 
deux dimensions r§cup£r£es par le systeme inf ormatique . 

25 La figure 2A repr^sente le cas d'une seule prise de vue 

dans un plan PI, fournissant une image II de l'os 1. La position 
de la source 7 a ete symbolis^e par un point oil convergent des 
rayons r de r6troprojection des quatre coins du plan PI. 

La figure 2B reprisente le cas d'une triple prise de 

3 0 vues dans des plans PI, P2 et P3, fournissant trois images II, 12 
et 13 en deux dimensions de l'os 1. En figure 2B, la position de 
la source n'est pas constante, celle-ci est diff&rente pour cha- 
que prise de vue. Toutefois, grace au localisateur 4, toutes les 
positions de la source sent connues dans le r^ferentiel de refe- 
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rence. Par soucis de clarte, seuls les rayons rl des quatre coins 
du plan Pi dans la position de 1 ! image II ont et§ repr6sent6s. 

On notera que la source de prises de vues peut subir, 
entre deux images en deux dimensions, d'autres mouvements que 
dans un meme plan comme illustre par les figures 1 et 2B. En 
d'autres termes, le bras 8 du systSme de prises de vues 6 peut 
presenter plus de deux degr6s de liberte, dent chacun pourra etre 
muni d'un dispositif de codage angulaire. Selon un mode de reali- 
sation prefere, le bras 8 est monte sur deux axes de rotations 
horizontal et vertical et un axe de translation vertical. 

On notera egalement que les images en deux dimensions 
peuvent n»£tre que partielles. Par exenple, dans 1 'application 
aux radiographics , celles-ci peuvent etre interpretees par I'ope- 
rateur pour valider les contours des vues en deux dimensions a 
prendre en compte pour la reconstitution. tfrie telle interpreta- 
tion n'est pas genante en raison du faible nombre de vues en deux 
dimensions n§cessaire selon l 1 invention (moins de dix en gene- 
ral) . La validation des contours dans le systeme informatique 
peut etre effectuee, par exenple, au moyen d'une souris, d'un 
crayon optique, d'un ecran tactile, ou equivalent, de fa<?on clas- 
sique pour un enregistrement de contour sur une image en niveaux 
de gris. 

Selon un autre mode de realisation pr6f6rS, la determi- 
nation des contours sur les images radiographiques est automati- 
see en mettant en oeuvre un procede dit de recalage en deux 
dimensions* Un tel proc§de consiste k determiner le contour de 
fa<?on autcmatique par analyse des gradients les plus forts en 
niveaux de gris. Ce procede est decrit, par exemple, dans la 
these de Gelu Ionescu, presentee publiquement le 4 decerribre 1998 
a l'Universite Joseph Fourier de Grenoble (France), et ayant pour 
titre "Segmentation et recalage d' images echographiques par uti- 
lisation de connaissances physiologiques et morphologiques ■ . 
Selon 1' invent ion, on met en oeuvre ce procede de recalage en 
deux dimensions en cocribinaison avec une projection, sur 1 ! image 
en deux dimensions, du module statist ique deformable de l'dbjet 
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dont on souhaite reconstituer les contours surfaciques en trois 
dimensions. On effectue une analyse iterative d'un contour pro- 
jet6 du modele sur 1' image en deux dimensions jusqu'a ribtenir, 
apres plusieurs passes (par exemple de 3 Sl 10) , une identity 
5 entre le contour du modele projet6 qui sera d6crit ci-aprds et 
dont 1 ' orientation et la d§formation auront ete adaptes par le 
procede de 1' invention, et le contour determine par analyse des 
gradients les plus forts en niveaux de gris. Cette combinaison 
sera mieux comprise par la suite en relation avec la figure 5. 

10 Pour simplifier la presente description, on consid^rera 

ci-apr£s que les contours des vues en deux dimensions sont 
connus. En pratique, les §tapes du procede de l 1 invention seront 
effectuSes en boucle, y ccmpris le recalage en deux dimensions, 
jusqu'3. obtenir un resultat satisf aisant . 

15 Selon l 1 invention, une fois que la ou les vues en deux 

dimensions ont et§ dbtenues, il reste a determiner, par recherche 
de correspond ance dans la base de donnfees, la forme et la taille 
de l'objet (ici, 1'os) en trois dimensions. 

Pour ce faire, l 1 invent ion pr^voit de rechercher le 

2 0 module pour lesquelles les distances separant les contours de 

chaque vue en deux dimensions et la surface du modele sont mini- 
males. 

Dans le mode de realisation prefSr6 oil on utilise un 
modele statistique d^formable (definissant ainsi une famille de 
25 moddles) de l»objet £ reconstituer, on recherche 3. s'approcher le 
plus possible, par dSformation iterative de ce modele, de la 
forme dont les contours sont tels que, corrnie on le verra par la 
suite en relation avec la figure 5, les rayons de retrqprojection 
s'appuyant sur les points des contours de 1' image (ou des images) 

3 0 en deux dimensions depuis la source (ou depuis leurs sources 

respectives) sont tous tangentiels £ la surface du modele, etant 
entendu que les differents contours en deux dimensions sont 
reperables dans le r^ferentiel de reference. 

La realisation du modele statistique deformable en lui 
35 m§me ne fait pas l'objet de la presente invention et est parfai- 
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tement classique. On se bornera a rappeler que la realisation 
d'un tel modele fait general ement appel a une recherche d'une 
forme moyenne d'une population d'obj ets de meme type (par exem- 
ple, des f&nurs) constituant la base de donn§es, suivie d'une 
5 analyse en composantes principales pour d6terminer les deforma- 
tions principales (essentielles) k appliquer a cette forme 
moyenne et cbtenir ainsi le tnodSle statistique. 

Une methode de determination d'un modele statistique §l 
partir d'une famille d'gchanti lions est d^crite, par exenple, 
10 dans 1 1 article "Building a Complete Surface Model from Sparse 
Data Using Statistical Shape Models: Application to Computer 
Assisted Knee Surgery" de Markus Fleute et St6phane Lavallee, 
paru dans MEDICAL IMAGE COMPUTING AND COMPUTER-ASSISTED 
intervention - MICCAI'98, Springer- Ver lag LNCS Series, pages 880- 
15 887, en Octobre 1998, et dont le contenu est incorporg ici par 
reference. On notera que la precision de la reconstruction des 
contours surfaciques qp£ree par 1' invention depend des e chant il- 
lons utilises pour construire la base de donnges du module sta- 
tistique. Par exenple, si ion moddle statistique de femur est cr6£ 

2 0 Sl partir d'une population d' §chantillons ayant des formes norma- 
les (sans pathologie) , on pourra recanstituer avec precision des 
formes normales rnais la precision sera limit^e si le f6mur radio- 
graphi6 a une forme pathologique . Par contre, si le moddle stati- 
stique est cr66 £ partir d'une grande population d 1 £chantillons 

25 contenant 3. la fois des formes normales et des formes pathologi- 
ques, on pourra reconstituer, avec precision, des dbjets normaux 
et pathologiques . 

Selon 1 ■ invention, une fois que l'on dispose des images 
en deux dimensions, on recherche le cas 6ch6ant & quel modele 

30 contenu dans la base de donn§es ces images se rapprochent le plus 
(par exemple, si plusieurs types d'os ou si plusieurs modeles 
statistiques d'un m§me os sont disponibles) , puis on applique des 
transformations k ce meddle pour s'approcher le plus possible 
d'une image en trois dimensions, telle que les rayons de retro- 

35 projection des images jusqu'3. leurs sources respect ives sont tan- 
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gentiels aux surfaces de l 1 image en trois dimensions, plac6e sur 
le trajet de ces rayons et dont ils sont respect iventent les plus 
proches . 

De preference, on commence par faire subir an module 
5 initial au repos, c ! est-a-dire correspondant a la forme moyenne, 
une transformation dite rigide, c'est-a-dire, ne touchant qu'a 
son orientation et sa translation dans l'espace. Puis, quand on a 
dbtenu la position pour laquelle les rayons de rgtroprojection 
sont tous ^ une distance minimale des contours du noddle initial, 

10 on fait subir une deformation non-rigide au modele, c'est-a-dire 
que sans modifier son orientation, on modifie sa forme 3. partir 
des donn^es contenues dans la base statistique en ttodifiant les 
coefficients des modes principaux du modele, jusqu'a obtenir la 
forme dont les contours se rapprochent le plus des rayons de 

15 r6troprojection. Le cas §ch§ant, le noddle est proposition^ de 
fagon tres approximative par l ! qp£rateur par rapport aux vues en 
deux dimensions af f ichees a l'ecran. On notera que les moyens 
(informatiques) de traitentent d 1 images utilises sont classiques 
dans leur structure et n ! ont done pas besoin d'etre d£taill§s. 

2 0 Les figures 3A a 3C illustrent ces deux Stapes du pro- 

cede de 1 'invention. La figure 3A represente une vue en trois 
dimensions d'un modele 21 avant toute deformation. II s'agit, par 
exenple, d'un modele statistique de 1 ' invention positionnS par le 
praticien dans une orientation approximative. La figure 3B repre- 

2 5 sente le modele 21 1 de la figure 3A a l f issue de l'6tape d' orien- 

tation sans modification de forme. La figure 3C represente 
1' image 22 en trois dimensions resultant de la mise en oeuvre de 
1' invention, c 1 est-a-dire correspondant au modele 21' , deforms 
pour que les rayons de r^troprojection soient (idSalement) tan- 

3 0 gents a tous ses contours. 

Par souci de sinplif ication, un seul groupe de rayons r 
de rStroprojection depuis une source 7 ont 6te repr§sent§s aux 
figures 3A a 3C et l'ecran de projection de l 1 image correspon- 
dante n»a pas ete represente. Comme il ressort des figures 3A a 
35 3C, un grand nombre de rayons r de la figure 3A ne sont pas tan- 
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gents aux contours du module 21. Le nombre de rayons de retrcpro- 
jection utilise depend de la pr6cision souhaitee et des caracte- 
ristiques du module statistique, essentielleraent, du nombre de 
points surfaciques choisis pour servir de reference lors de la 
definition de ce module statistique. 

De la figure 3A a la figure 3B, la transformation est 
"rigide", c'est-a-dire que le modele n'est pas d§f ormS . Lors de 
cette transformation, le module 21 subit des translations et 
rotations afin d'obtenir des rnesures d' hearts minimum entre 
chaque rayon de retroprojection partant des points de contour de 
I 1 image et la surface de l'objet dans sa position courante . On 
utilise pour cela des methodes mathematiques de recherche du 
minimum de la somme des carres de ces ecarts par rapport aux six 
paramStres (trois degres de liberte en rotation et trois degr£s 
de liberty en translation) d§finissant la transformation rigide 
recherch6e. On obtient une configuration telle qu»illustree par 
la figure 3B dans laquelle le module 21 1 a une position et 
orientation correcte mais une forme encore imparf aite . C'est 
pourquoi certain rayons r' ne sont pas tangents mais t raver sent 
le modele. 

De preference, la position resultant de la transforma- 
tion rigide (rotation, translation) est obtenue par une mSthode 
(algorithme) dite du point le plus proche (Iterative Closest 
Point, ICP) . Un exemple d'une telle methode est d^crit dans 1' ar- 
ticle "A Method for Registration of 3-D Shapes" de P au l J. Best 
et Neil D. McKay, paru dans IEEE TRANSACTIONS ON PATTERN ANALYSIS 
AND MACHINE INTELLIGENCE, vol. 14, NO. 2, en fevrier 1992, et 
dont le contenu est incorpore ici par reference. 

De la figure 3B a la figure 3C, la transformation est 
"non- rigide" (eiastique) et consiste Sl deformer le module 21' 
sans changer son orientation jusqu'& obtenir 1 1 image 22 en trois 
dimensions pour laquelle tous les rayons de retroprojection r" 
sont a une distance minimale des contours (idealement tangents) . 
Cette transformation non-rigide est effectu6e, par exemple et 
selon un premier mode de realisation, en calculant les paramStres 
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de deformation a partir d'un algorithme connu sous la denomina- 
tion anglo- saxonne "Down Hill Sinplex" . Cet algorithme a §te 
d§crit, par exemple, en 1965 par J. A. Nelder et R. Mead dans 
COMPUTER JOURNAL, vol. 7, pages 308-313, dont le contenu est 
5 incorpore ici par reference. 

Selon ion mode de realisation pre fere de 1 1 invention, 
l 1 algorithme utilis§ pour determiner la transformation non-rigide 
du modele statistique est basS sur la methode de Levenberg- 
Marquardt. Le principe de cette m£thode est de determiner le 

10 minimum d'une fonction non-lin§aire multidimensionnelle en utili- 
sant les derivees part ie lies de la fonction par rapport aux para- 
metres de deformation du modele. Cette methode a £t6 decrite, par 
exemple, en 1963 par D.w. Marquardt dans JOURNAL OF THE SOCIETY 
FOR INDUSTRIAL AND APPLIED MATHEMATICS , vol. 11, pages 431-441, 

15 dont le contenu est incorpore ici par reference . La fonction dont 
on cherche le minimum est ici la scmme des carr6s des distances 
entre un ensemble de rayons de r6troprojection s'appuyant sur les 
points de contours et la surface du modele. Les parametres de la 
recherche du minimum sont les coefficients que l ! on applique a 

20 chaque mode de deformation du modele statistique, ainsi que cela 
est pr6sent£ dans l 1 article "Building a Complete Surface Model 
from Sparse Data Using Statistical Shape Models: Application to 
Conputer Assisted Knee Surgery" de Markus Fleute et St£phane 
Lavallee, d§ja cite. 

25 Cette methode est appliqute et repetee a chaque 

position approchee jusqu'a ce qu'il ne soit plus possible de 
minimi ser les hearts. On enchaine ainsi un recalage rigide, un 
recalage elastique, un recalage rigide, un recalage elastique, 
etc. Une fois la forme en trois dimensions reconstitute, sa 

3 0 taille est bien entendu egalement connue dans la mesure ou toutes 
les vues en deux dimensions ont des tailles connues, de meme que 
la position, dans le referent iel de r6f6rence, du module par rap- 
port aux sources. Le nombre de passes necessaire depend essen- 
tiellement du nombre de parametres caracteristiques utilises pour 

35 mod§liser la deformation. 
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On notera que l»ordre dans lequel sont effectuees les 
deux etapes ci-dessus est inportant. En effet, si on commence par 
effectuer une deformation non-rigide, on va faire correspondre un 
modele avec une mauvaise orientation qu'il sera tr§s difficile a 
r6cup£rer en raison de la deformation de forme subie. Toutefois, 
la minimisation des hearts peut etre effectu^e gldbalement en 
prenant en cortpte les paramStres rigides et les param&tres de 
deformations en meme temps §l chaque iteration du processus mat he - 
tnatique de minimisation. 

On notera que, k chaque etape iterative du proc§d§ dans 
laquelle on calcule les hearts entre les contours du modele et 
les rayons de retrqprojection, de nouveaux points entrant dans la 
raesure sont determines cote modele et/ou cote rayons de projec- 
tion. Ceci constitue une distinction import ante par rapport & des 
proc£d6s connus de reconstitution d 1 images en trois dimensions 
dans lesquels les points de mesure sont les m§mes pour toutes les 
iterations . 

Utie autre caracteristique de la presente invention est 
que les evaluations des ecarts minimaux ne sont pas necessaire- 
ment faites de fagon exhaustive pour tous les elements de surface 
du module. En effet , cela conduirait parfois 3. des calculs trop 
dispendieux meme si le nombre d ! elements de surface peut etre 
minimise en utilisant des elements triangulaires . Selon 1 ' inven- 
tion, les evaluations d 1 ecarts sont faites par rapport a des 
lignes caracteristiques constituant des generateurs de contours 
du moddle. Cela minimise le nctribre de points de mesure. Par exam- 
ple, un os de type femur, defini par environ 5000 bords de surfa- 
ces triangulaires, est defini par environ 300 generateurs de 
contours . 

De preference, 1» analyse effectuee par le procede de 
1» invention ne touche que des generateurs de contour de l 1 image 
en trois dimensions. Ainsi, la surface du modele statistique est, 
de preference, formee d 1 elements triangulaires dont certaines 
arites definissent des generateurs du contour. Un generateur de 
contour est defini par les aretes (bords) des triangles qui, pro- 
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jetes sur un plan, definissent le contour (interne ou externe) . 
Le recours aux g6n6rateurs de contours permet de reduire conside- 
rablement (par exemple, d'au moins un facteur 10) le ncmbre de 
points a chercher dans le modele pour verifier la correspondance 
5 avec 1' image reconstitute. 

La figure 4 illustre la definition d'un g6n6rateur de 
contour d'une image (meddle) en trois dimensions dont la surface 
est formee d' elements de surface triangulares tels que dgcrits, 
par exemple, dans 1' article "Anatomy-based registration of ct- 

10 scan and intraoperative x-ray images for guiding a surgical 
robot" de A. Gueziec, paru dans IEEE TRANSACTIONS ON MEDICAL 
IMAGING, 17(5), pages 715-728, en octofore 1998, et dent le 
cOTitenu est incorpor£ ici par reference. 

La figure 4 reprSsente, de fa<?on siraplifi6e, deux 

15 triangles 30 et 31 def inissant une portion de surface d'un modele 
en trois dimensions et dont 1' arete commune constitue une g6nera- 
trice 32 du contour. Le calcul a effectuer pour determiner si une 
arSte constitue ou non une g#n£ratrice du contour consiste cl cal- 
culer les angles a et p respectifs entre les normales aux surfa- 

2 0 ces des triangles 30 et 31 et les rayons 33 et 34 reliant le 

centre de ces surfaces au centre 35 de projection correspondant a 
la position de la source de rayons X. Si l'un des angles est 
inferieur a 90° tandis que 1' autre est sup6rieur a 90° , leur 
arete cotnnune 32 est alors une g§n£ratrice du contour. 
25 La figure 5 illustre le type de mesure effectu^e pour 

determiner l'ecart e entre un rayon r de r^troproj ection et un 
point d'un generateur de contour. Cette figure represente, de 
fa<?on sch£matique, une image I de 1'bbjet dans un plan P et une 
forme 21 du modele statistique placSe sur le trajet des rayons de 

3 0 r^trqproj ection de 1» image I jusqu'd la source (non representee) . 

Pour chaque rayon de retrcproj ection r, on cherche les 
points respectifs 41 et 40 du rayon de retroproj ection et d'un 
generateur de contour choisi parmi tous les generateurs, pour 
lequel l'Scart e entre ces points est minimal. 
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Un avantage de la prgsente invention est qu'elle permet 
une reconstitution beaucoup plus rapide d'une image en trois 
dimensions par rapport aux techniques connues. 

Un autre avantage de la pr£sente invention est qu'elle 
5 permet un alignement correct du modele mime dans des zones oQ on 
ne dispose pas de donnges de projection par les images en deux 
dimensions (par exemple, certaines courbures internes) . 

On notera que le contour de surface pourra, le cas 
£ch6ant, §tre affin£ au moyen d'un palpeur m§canique, qptique ou 
10 magnet ique fournissant des coordonn£es spatiales dans le m£me 
r£f§rentiel que les images en deux dimensions. On obtient alors 
des points supplementaires qui peuvent §tre utilises dans la 
recherche de 1» image en trois dimensions. Le recours 3. un palpeur 
peut servir, par exerrple, & diminuer le nombre de vues en deux 
15 dimensions ngcessaires en donnant une information en trois dimen- 
sions. 

On notera egalement que le procedS de l f invention peut 
s'appliquer & plusieurs contours surfaciques en trois dimensions 
constituant un ou plusieurs objets. Par exertple, on peut ainsi 
2 0 recanstituer sitnultanement la surface externe et la surface cor- 
ticale interne d'un os. On peut ggalement reconstituer la forme 
de plusieurs os participant a un articulation et visibles simul- 
tanement sur les images radiologiques (par exemple, le tibia et 
le f§mur si on s'int§resse au genou) , ou encore la forme de plu- 

2 5 sieurs fragments osseux d'un raeme os. Le module recherchS doit, 

dans ce dernier cas, contenir les transformations rigides entre 
chaque contour surfacique en trois dimensions. Dans d' autre cas 
d' application Sl plusieurs objets, il faut bien entendu connaitre 
les differents modules statistiques et les gventuelle transforma- 

3 0 tions rigides ou 6lastiques entre eux. 

Bien entendu, la pr£sente invention est susceptible de 
diverses variantes et modifications qui apparaltront k l'homme de 
l'art. En particulier, le nombre de vues en deux dimensions 
devant §tre utilise pour la reconstitution depend de la precision 
35 souhaitSe et de la complexity du modele. Dans certains cas, une 
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seule vue en deux dimensions pourra suffire. De plus, bien que 
1 ! invention ait ete decrite en relation avec une source de rayons 
X mobile, celle-ci pourra €tre remplac£e par plusieurs sources 
fixes, pourvu de respecter la contrainte d'obtenir des images en 
5 deux dimensions dans un mime ref erentiel . En outre, la mise en 
oeuvre de 1' invention fait bien entendu appel 4 des techniques de 
traitement num£rique d 1 images dont la r§alisation pratique est a 
la portee de l'homme du metier a partir des indications fonctian- 
nelles donnSes ci-dessus et dans les publications cities en ref6- 
10 rences. 
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REVENDIC&TIONS 

1. Proc£de de reconstitution d'une image en trois 
dimensions repr£sentant les contours surfaciques d'au moins un 
objet (1), k partir d'au moins une vue en deux dimensions de cet 
cibjet prise aux rayons X, caract6ris£ en ce qu'il consiste Sl : 
5 determiner la position de la source (7) de prise de 

vues dans un rSfgrentiel de r6f£rence ; 

select ionner un modele predefini const ituant une forme 
moyenne (21) de 1' objet / et 

de f agon iterative jusqu ■ £ ce que les contours du 
10 modele soient tels que les ecarts, entre des rayons de retropro- 
jection des contours de 1 ! image en deux dimensions depuis la 
source et la surface du modele , soient minimaux : 

selectionner une orientation et une position du meddle 
dans le r£f£rentiel de reference, puis 
15 selectionner vine deformation du modele pour modifier 

ses contours en trois dimensions. 

2. Proc6d§ selon la revendication 1, caracterisS en ce 
que le moddle est obtenu 3. partir d'une population d'objets pour 
laquelle on recherche la correspondance statistique commune §l 

2 0 tous les cbjets pour determiner une forme moyenne et les deforma- 

tions principales par rapport & cette forme moyenne, de fagon £ 
disposer d'au moins un modele statistique. 

3. Proc6d6 selon la revendication 2, caracteris£ en ce 
que les etapes de selection iteratives consistent Si faire subir 

25 au modele statistique, successivement, une transformation rigide 
modifiant sa position et/ou son orientation et une deformation 
non- rigide modifiant ses contours surfaciques. 

4. Proc£d§ selon l'une quelconque des revindications 1 
a 3, caracterise en ce que les contours de 1' image en deux dimen- 

3 0 sions sont obtenus automatiquement en projetant le modele dans le 

plan de 1' image en deux dimensions, et en d£formant les contours 
projet^s de fagon d les faire coincider avec les points de fort 
gradient en niveaux de gris de 1' image en deux dimensions. 
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5. Procede selon la revendication 4, caracterise en ce 
que la determination automatique des contours de 1' image en deux 
dimensions est effectu6e de fagon iterative, chaque iteration 
etant intercalee entre deux iterations successives des Stapes de 

5 selection. 

6. Procede selon l'une quelconque des revendi cations 1 
a 5 r caracterise en ce qu'il consiste & determiner, dans le rgf§- 
rentiel de reference, des coordonnees en trois dimensions de 
points de l'objet, de fagon k disposer de points de reference 

10 supplementaires pour les etapes de selection iteratives de posi- 
tion, orientation, et deformation. 

7. Procede selon l , une quelconque des revendi cations 1 
Sl 6, caracterise en ce qu ! il consiste & utiliser plusieurs images 
en deux dimensions pour lesquelles les positions respectives de 

15 la source de prise de vues sont toutes determin^es dans le r£f6- 
rentiel de reference, et §. effectuer les Stapes de selections 
iteratives en tenant compte des rayons de retrqprojection des 
contours de toutes les images en deux dimensions. 

8. Procede selon la revendication 7, caracterise en ce 
20 que le nombre d 1 images utilise est fonction de la precision 

souhaitee . 

9. Procede selon I'une quelconque des revendications 1 
a 8, caracterise en ce que la surface du modele est constituee 
d 1 elements de triangle, lesdits e carts etant mesures par rapport 

25 a des points de certaines aretes constituant des generateurs du 
contour en trois dimensions. 

10. Procede selon l'une quelconque des revendications 1 
i. 9, caracterise en ce qu'il est applique k la reconstitution des 
contours surfaciques de plusieurs objets lies entre eux par des 

3 0 relations de transformation rigide et/ou eiastique. 

11. Procede selon l'une quelconque des revendications 1 
a 10, caracterise en ce qu'il est applique k la reconstitution 
d' images d'os. 
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12. Systdme de traiteraent d' images, caract§ris€ en ce 
qu'il ccxtporte des moyens pour la mise en oeuvre du proced6 selon 
l'une quelconque des revendi cations 1 k 11. 
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